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We review some of the well-known features of quantum cosmology, such as the factor ordering 
problem, the wave function and the density matrix, for a dark energy dominated universe, where 
analytical solutions can be obtained. For the particular case of the phantom universe, we suggest a 
quantum system in which the usual notion of locality (non-locality) of quantum information theory 
have to be extended. In that case, we deal also with a quantum description where the existence of 
a non-chronal region around the big rip singularity is explicitly accounted for. 
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I. INTRODUCTION 

Besides confirming that the universe expands in an 
accelerated fashion, recent data coming from SNIa and 
other observations 1][2] leave the issue of the precise way 
in which such an accelerated behavior actually occurs un- 
settled. In fact, the possibility for a superaccelerated ex- 
pansion beyond what is predicted by a cosmological con- 
stant has been raised, implying serious theoretical dif- 
ficulties. Several models have been proposed so far in 
order to explain accelerated expansion. The most pop- 
ular among them are the so-called quintessence models, 
which are characterized by a universe filled with an ho- 
mogeneous fluid with an equation of state p — wp, where 
p and p are the pressure and the energy density of the 
fluid, respectively, and w is a constant within the range 
— 1 < w < — |, being the value of w = — 1 equivalent to 
a cosmological constant. This would trigger by itself an 
exponential expansion if one assumes no other energy or 
matter sources in the universe. This acceleration would 
correspond to a value for w less than — 1 , so allowing for 
what is called the phantom regime. This regime entails 
violation of the dominant energy condition, and might 
imply some interesting features from the point of view of 
the quantum theory. Some of such features will be an- 
alyzed in this paper within the realm of a more general 
formalism where other possible scenarios are also con- 
templated. 

On the other hand, from a quantum mechanical stand- 
point, the universe is a rather special system since it 
cannot be described as a whole in terms of space-time 
coordinates but in terms of geometries. Hence, it offers 
a particularly interesting framework to deal with some 
well-known issues of quantum mechanics, such as those 
related with the notion of non-locality. 

In this paper we shall consider therefore the quantum 
theory of an accelerating universe which is filled with 
dark energy, both when the dominant energy condition 
is satisfied and for vacuum contents where such a condi- 
tion is manifestly violated. In the latter case the notion of 
non-locality is discussed in a multiverse scenario, where 
it must be necessarily generalized or extended. The gen- 
eralized quantum theory of Hartle Q is then applied to 



cases where the future singularities are replaced for non- 
chronal bounded regions. 

The paper can be outlined as follows. Sec. II deals 
with a phantom universe covering the entire time inter- 
val in such a way that it becomes describable as a mul- 
tiverse. In sec. Ill we review the canonical Hamiltonian 
formalism for a quantum universe, particularizing in the 
problems related with the density matrix and the possi- 
bility for the existence of entangled states in the phantom 
multiverse. The generalized quantum theory is applied 
to the quantum multiverse in which the big rip singu- 
larity is replaced for a bounded non-chronal region in 
sec. IV, where a decoherence function is used and ob- 
servable probabilities are isolated from it using the Har- 
tle procedure. We summarize and conclude in sec. V. 
An appendix on the orthogonality properties of the wave 
function is also added. 



II. THE PHANTOM MULTIVERSE 

For a flat Friedmann-Lemaitre-Robertson- Walker 
(FLRW) universe filled only with dark energy, the equa- 
tions of motion can be obtained from the Hamiltonian 
constraint in the phase space, 
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where p a the canonical momentum, G is the gravitational 
constant, po is the energy density at the coincidence time, 
i.e., the time in which the dark energy started to dom- 
inate the expansion of the universe, and a is measured 
in units of ao- In deriving Eq. ([1]) we have used the 
integrated form of the energy conservation law, 
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In the configuration space, the corresponding Friedmann 
equation reads 
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and hence the scale factor runs as 

a(t) oc (t br ± t)~3(|J|-i) 
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where ifc r is a constant, and the + and — signs stand for 
the quintessence and phantom regimes, respectively. The 
energy density goes then as 
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so, in the phantom case, t\, r turns out to be the time at 
the so-called big rip singularity [|[ , where the scale factor 
and the energy density blow up to infinity (see Fig. [I]). 

If one had to look at this phantom universe as evolving 
from an initial coincidence time to infinity, then in order 
for the scale factor, a, to be positive in (fj|, not all values 
of the parameter w entering the equation of state for 
the universe would be allowed. Actually, for a phantom 
regime only the values given by [3| 
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are allowed as the scale factor becomes ill-defined at 
t > tt, r otherwise. Moreover, for observers staying at 
times before the big rip the singular character of this 
would necessarily imply that the big rip should be cut off 
from the considered physical manifold, so making com- 
pletely unphysical the region beyond the big rip singu- 
larity. Thus, besides introducing the condition (JSJ) one 
must also take into account either that stable wormholcs 
and ringholes (whose existence is induced by the violation 
of the dominant energy condition implied by the phan- 
tom nature of the cosmic fluid) would crop up at both 
sides of the singularity shortcutting the space-time Q , so 
keeping the singularity outside the trajectories followed 
by physical signalling and making therefore accessible to 
any observers the contracting region beyond the big rip, 
no matter whether it is cut off or not; or even by con- 
sidering that quantum effects could somehow smooth out 
the singular character of the big rip. The condition © 
could be thus interpreted as a condition for the existence 
of a multiversal scenario in which every value of n would 
provide us with a different universe with its own time 
evolution for the scale factor. It is worth considering 
that this multiversal scenario can still be a classical one, 
though its quantum counterpart presents some interest- 
ing features that we will discuss later on. 



III. THE QUANTUM DARK ENERGY 
UNIVERSE 

A. The wave function 

Following the canonical quantization procedure, the 
dark energy dominated universe would be quantum- 
mechanically characterized, in the configuration space, 
by a wave function, $(a), which is annihilated by the 




FIG. 1: Wormholes could connect both sides of the universe 
avoiding that way the big rip singularity. For times t,ti < t < 
tf, i.e., the achronal region, the manifold cannot be foliated. 



Hamiltonian density fl}, m its operator form, satisfying 
in this way the quantum version of the Hamiltonian con- 
straint. But in doing that we have to make a particular 
choice of the factor ordering between the conjugate vari- 
ables. There are many ways in which one can depict this 
ambiguity, maybe a quite general one would be the fol- 
lowing. Classically, the conjugate variables commute so 
the rhs and lhs in the following expression are equivalent 
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but, of course, their quantum counterparts are not nec- 
essarily equivalent as they depend on the commutation 
relations between the conjugate variables. For instance, 
for the kinetic term of the Hamiltonian (TTJ) , we would 
obtain 
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and a and j3 are related to the exponents in Eq. ([7]) 
through 
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so under canonical quantization, p a — > —ihJ^, a and 
(3 represent the factor ordering ambiguity in the cor- 
responding Wheeler-DeWitt equation |6j, which for the 
case being considered is 
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where A 2 , = -^qPo, and N is the lapse function. Here, as 
well as in Eq. }T]), we have assumed: (i) that dark energy 
dominates in such a way that all matter in the universe is 
subdominant and can therefore be disregarded, and (ii) 
that the quantum state is given as a pure state. The 
first assumption stems from the feature that the 70% 
of the universal energy is made up of dark energy, even 
though other fields should be considered in a more de- 
tailed framework. The second assumption can be thought 
of as a plausible one in spite of the state derived from it 
is not the most general. The solutions of the Wheeler- 
DeWitt equation, Eq. (fT3i|) , can be given in terms of 
Bessel functions 

*(o) = a a C„(\a q ), (13) 

where 

q =l(l- w ),\=* v\^=c?-(5\ (14) 
2 nq 

C is the Bessel function of the first or second kind, J v 
and y v , respectively, and w satisfies the discretization 
([5]) for the case of the phantom multiverse, and where 
we have kept the Planck constant. This wave function 
would quantum-mechanically describe a quintessence en- 
ergy dominated universe, and actually the state of the 
quantum universe for the phantom regime, were it not 
for the possible existence of a noncausal multiply con- 
nected region around the big rip singularity in the phan- 
tom energy dominated universe, at least for the fraction 
of physical reality which always be outside the connec- 
tion. That wave function must satisfy, as usual, given 
boundary conditions. For these boundary conditions we 
choose: (i) the wave function ought to be regular every- 
where, even when the metric degenerates at time t — ► oo 
in the phantom universe, and (ii) it should vanish at the 
big rip singularity when a — > oo. These conditions are 
satisfied by the wave function (|13p if we enforce a to be 

<*<§, (15) 

for the quintessence regime for w > — 1 and for w < — 1 
just in the region before the big rip, t^.. In the latter case, 
for times after the singularity we must add the condition 

a±qRe{v)>0, (16) 

where the + sign stands for the Bessel function of the first 
kind, J v , and the — sign does for the Bessel function 
of the second kind, y v . A linear combination of these 
two solutions together with the boundary conditions (|15[) 
and (|16p would represent, therefore, the quantum state 
of a dark energy dominated universe. From this state we 
should be able to recover the semi-classical universe in 
which we live. This may be accomplished by taking the 
limit h — > in the expression for the wave function. In 
particular, using the asymptotic expansion of the Bessel 
functions [7J for the wave function (|13p . we can obtain 



that in the semiclassical approximation 

/_^_ a ("-29) e ± ''( Aa ' ! -5 !y7r -3 7r ). (17) 
V 7tA 

This wave function represents the state of the classical 
universe in the sense that it is a quasi-oscillatory wave 
function whose argument is essentially the classical action 
(So = ^ a9 li so that the correlations between the classical 
variables [g] are satisfied, i.e, p a — where p a is the 
classical momentum, is the equation of motion, and 

A = a a ~i q (18) 

in Eq. (fTT|) is a prefactor smooth enough to satisfy the 
Hartle criterion [9J. In fact, 

^t|-(|A| 2 ^)=0, (19) 

which in our case implies 

- #- (a 2a ~^) ~ a 2 "" 3 - 0, (20) 
a oa \ oa J 

that is, the Hartle criterion is satisfied for quintessence 
models (w > —1) for the boundary condition (fT5")) . For 
the phantom regime we should replace 2a < q for 2a < 3 
for the boundary condition, so that Eqs. (fT5|) and (|2"U|) 
are both satisfied. In both cases, the semiclassical ap- 
proximation (|17p obeys the Hamilton- Jacobi equation in 
the limit h — > 0, irrespective of the value of q in the 
phantom multiverse that is described by ([6]), and for all 
choices of factor ordering satisfying the boundary condi- 
tions. In particular, any universe in this phantom mul- 
tiverse would have a semi-classical domain described by 

(HZ}. 

B. The density matrix 

The most general quantum state would be given how- 
ever in terms of a mixed density matrix rather than a 
pure wave function [lfj. One can also compute the den- 
sity matrix for the case being considered by taking 

/>oo 

p(a , ,a)= dTK(a',T;a,0), (21) 
Jo 

where K(a! , T; a, 0) is the Schrodinger propagator and 
the integration over time is introduced to account for 
the invariance of the time separation between any two 
hypersurfaces. 

In the case of our dark energy universe, we can take the 
gauge J\f = a 3 in Eq. ([T2"]) , so we get a set of Hamiltonian 
eigenf unctions 

iT$ fc (a) = Pi $ fc (o) (22) 

given by 

S fc (a) = N k a a J k {\a"), (23) 
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where Nk is a normalization factor, and q and A are given 
by Eq. (fT4")l . with eigenvalues 
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weighted by the function W(a) = a~' 2a+1 \ the set of 
eigenfunctions would satisfy the following normalization 
relations 



($fc(a)|$ fc (a)) = l,V*>0, 



in which we have used in Eq. ([23)1 Nk 

way, for k ^ I, we obtain 



(26) 
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($fc(o)|$i(o)) = ,{k-l) even 

($ fc (a)l^(a)) = | ( " 1) ^:7 )Ve ,(fc-0 odd. 

This set of equations can be considered as orthogonality 
relations in the sense that they permit to split the whole 
Hilbert space spanned by the Hamiltonian eigenfunctions 
into two Hilbert subspaces, i.e., the subspaces spanned by 
the odd and even modes. In that case, the Hamiltonian 
eigenfunctions form two orthogonal basis, in the usual 
sense, for the subspaces. 

But the zero mode is not normalizable. We may reg- 
ularize it by using some cut-off or minimum length, l p , 
taking on the limit l p — > at the end of the calculations. 
We have (see the Appendix), 
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with which we could take the normalization relations 
for all k > 0, with a suitable normalization factor, Nq. 

With that regularization, we can make use of the fol- 
lowing set of functions 
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which are constructed in terms of the Laguerre polyno- 
mials, L n (x), 
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We have in this way obtained an orthonormal set un- 
der the scalar product (|2"5|) , which can be used as the 
basis for the square integrable functions upon which the 
Hamiltonian would act, so 



|*»)<*n| =Id. 



Using this set, we have 
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where eg = a 2 — (3 2 > for the parameters r and s in where 
Eq. CZJ to be real (see Eqs. (HO]) and (fTTji ). 

However, this set of eigenfunctions are not orthogonal. 
For instance, choosing for the scalar product 



D k i = y^CrefcC 1 *;, 



(33) 



and the Cy's are the coefficients for the change of basis, 
i.e., 

^n{a) = Y,C nm <5>m(a). (34) 

m 

Hence, the propagator can be written as 

K(a',T;a,0) = ^ D u {a'\e* TH \$ k )($i\a) 

kl 

= ^£> w e* T/3 ^ fc (a')^(a), (35) 
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and the density matrix computed to be 

$fc(o')^i(o) 

kl 



p{a,a) = ^D kl 

q t e 



(36) 



where, in order to make the integral well-defined, we have 
Wick rotated time counterclockwise. Wick rotating in 
the opposite direction would have implied inserting the 
identity before the evolution operator and setting a minus 
sign in the exponent. 

In the case of the dark energy universe, with the Hamil- 
tonian eigenfunctions (f2"3"]l . the coefficients in Eqs. (|3"3"|) 
and (j34|) can be computed. For, we can take advantage 
of the properties of the first kind Bessel functions, which 
form up an ovcrcompletc set in the sense that any ar- 
bitrary function can be decomposed in terms of Bessel 
functions through a Neumann's expansion [ll[ 



Z V fk{z) = ^2 CknJn+»(z), 
n=0 

with the coefficients being given by 

Ckn = -^—.( dt f k (t)A, hU (t), 

*™ J\t\<R 



(37) 



(38) 



where R is the distance from t = to the closest pole of 
fk(t), and the A ni „, the Gegenbauer's polynomials, are 
defined by 



m— x 

Thus, we can rearrange Eq. (|37|) so that any arbitrary 
function could be written as 



(31) g k (a)=a a (\a^ fk(\a q ) = J2 c kna a J n+ „(^ q ). (40) 
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In particular, with v = 5, the orthonormal set of func- stant, are given by Eq. (|38|) . with 
tions (|2T)j) can be decomposed as 



^k(a) = ^C fc „$„ + i(a) 



(41) 
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in which the coefficients, unless for a normalization con- that is, 
I 
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where m is a non-negative integer. The density matrix 
(l36l) can then be written as 



p(a',a) 



^ $ fc+ i(a')^+i(a) 



(44) 



A density matrix for a physical system is supposed to 
be definite positive, and Eq. (|44|) is however not neces- 
sary so, even for positive values of the coefficients Dki- 
Parameters a and [3 have nevertheless no clear physical 
meaning since a semiclassical state should be indepen- 
dent of the particular choice of a and (3. We then could 
still take for the physical state of the universe the re- 
duced density matrix resulting from integrating out a 
and j3, i.e., 

p r (a',a) = J da J d(3 p(a! ,a;a, (3), (45) 

This integral turns out to be divergent so that one could 
not obtain a meaningful density matrix for the state of 
the universe. One way to avoid this problem could be 
taking some particular values for the factor ordering as a 
boundary condition. For the particular choice a = (3 = 
we would in fact have 



p(a 
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(46) 



where the coefficients are given by Eq. (|33|) . However, 
even though this density matrix does not show the usual 
divergences due to vanishing values of the denominator, 
the coefficients D nm in it are still divergent. 

On the other hand, our particular choice leading to 
the density matrix (|46|) corresponds to a rather arbitrary 
choice of the factor ordering. Other choices could also 
imply a definite positive density matrix. It would fol- 
low that an alternate philosophy could be constructing a 



propagator in terms of pairs of levels instead of the single 
levels which correspond to the Hamiltonian eigenvalues, 
that is 



(a'|e iTff| 
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In that case, we have 

P(a , a = > D nm ^-^2 ^ 48 

■* — ' g z (m z — n z + m — n) 

nra v ' 

i.e., although the factor ordering ambiguity is no longer 
present in the denominator, the diagonal elements be- 
come now divergent. 

The kinds of divergences and unphysical states that we 
have just uncovered were already pointed out by Hawk- 
ing and Page for density matrix in quantum cosmology 
[l(j ] [e|. Actually, these difficulties can be seen to arise 
because the system can reach a state with zero value for 
its Hamiltonian eigenvalue, and would be expected to be 
solved in the framework of a proper quantum theory of 
gravity, in which a minimum energy, the Planck mass, 
m p , should exist. In such a case, a nonzero minimum 
value of the Hamiltonian eigenvalue would be expected 
that rendered the density matrix definite positive and 
always convergent. 



C. Entangled states in the multiverse 

Let us, now, be concerned with a phantom universe 
and its big rip singularity, then we could consider the 
case in which no wormholes are connecting the regions 
before and after the big rip; i.e. when the wormholes 
that branch off in the neighborhood of that singularity 
simply connect two asymptotic regions on the same side 
of the singularity. If thereby such wormholes are disre- 
garded and the hypersurface at the singularity is cut out 
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so that the whole space-time is divided into two separate 
parts, then there will be two independent wave functions 
which should be associated with different realizations of 
the boundary conditions and distinct time intervals. The 
first of these intervals runs from the coincidence time un- 
til the time at the big rip, and the second one goes from 
the latter time until infinity. The general boundary con- 
ditions that the quantum state be regular everywhere 
and exactly vanishes at the big rip singularity amount 
to a wave function which should be generally expressed 
in terms of a different linear combinations of first and 
second kind Bessel's functions, J and y, on each inter- 
val. These two wave functions for both sectors can be 
regarded to play the role of some bases for the quantum 
state of a specific n-phantom universe. So, in general we 
can describe this state as 

*„ = c n ^ r i + c'Ij^j, (49) 

with 

(c7f + (fii? = 1. (50) 

Let us interpret for a moment the integer number n 
defined in Eq. )6j as a quantum number labeling the 
different universes in our multiverse, and then consider 
the quantum states for two universes with different values 
of the quantum number n. In that case, ipj and ipj'j, 
which quantum-mechanically describe the regions before 
and after the singularity for a single universe labeled n, 
are both strongly peaked at time-like separated regions 
and can therefore be correlated to the similar regions of 
the other universe, say m, as all of such states satisfy 
the so called Hartle criterion ([20| . The result for the two 
universes could be a common state 

$ = c?cfV}V? + chcTi^nWi- (51) 
I 



D(a',a) = NTr 



where the P^'s are projection operators which forms up a 
set, {P ap }, that corresponds to the exhaustive and exclu- 
sive set of alternatives defined on a given non- intersecting 
spacelike surface E, and the Oj's are particular sequences 
of coarse-grained alternatives {a} = ai,...,otj that de- 
scribe particular histories. The exhaustive set of histories 
consists then all possible sequences {a}. The function 
p denotes the density matrix encompassing the bound- 
ary condition of the system on an initial nonintersecting 
spacelike surface So- It will be given either by the fac- 
torizable probability p = W = $(a)$(a') or by the ex- 
pressions for the density matrix described in the previous 
section, if the probability function W is not factorizable, 



Since the singularity has been cut off, mixed states tpftpfj 
and ip^jipj 1 are no longer possible if there are correla- 
tions between tpfipj 1 and ipfjip^i an d' therefore, state 
([51) should be an entangled state. 

Eq. (|5ip can straightforwardly be generalized to the 
infinite possible number of universes and would imply 
that knowing the state of our universe we would auto- 
matically know the state of the other universes belong- 
ing to the same multiverse scenario. This is a cosmic 
translation from what is usually dubbed quantum non- 
locality in quantum information theory. Since there is no 
space-time between any two universes in this multiverse, 
the term locality (or non-locality) used to characterize 
correlations between two particles in a common space- 
time, is no longer suitable to physically characterize the 
above correlations between universes, provided that lo- 
cality refers to just space-like or time-like location in a 
common space-time. It would instead refer to correla- 
tions between the quantum states of different universes, 
which would become entangled as a result. 

IV. GENERALIZED PHANTOM UNIVERSE 

However, the existence of a non-chronal region around 
the big rip makes it impossible to have any quantum 
state for the phantom universe. Actually, in order for 
having a proper quantum theory of one of the universes 
of the multiverse we have to make use of a generalized 
quantum theory 0. Technically, the physical system we 
ought to deal with consists of a space-time manifold con- 
taining an intermediate bounded non-foliating region on 
the neighborhood of the big rip singularity, filled with 
closed time-like curves (CTCs), which is chronologically 
placed between two regions that are both foliable by a 
family of nonintersecting spacelike surfaces E p (See Fig. 
1). We then introduce the generalized decoherence func- 
tion of Hamiltonian mechanics which reads 



(52) 

I 

i.e., for a mixed state. A is a generalized evolution ma- 
trix that can be defined in terms of a nonunitary matrix 
X s which replaces the usual unitary evolution matrix U 
of the decoherence function for fully foliable manifolds. 
It can be given by 

A = (7(E / ,E oo )- 1 A s [/(E l ,I] ), (53) 

in which Eo and Eoo are surfaces at the furthest possi- 
ble past and future, the latter being assumed to be at 
the event horizon (i.e. in the present case at infinity), 
and Ej and Ej are the latest and earliest nonintersect- 
ing spacelike surfaces, after and before the nonfoliable 
region, respectively. Finally, the normalizing factor N is 



P*,(X P )...P*£ ^ k+l )XP*{Y> k )...Pl,{i: l )pPl i ^^ 

P fc + 1 IB X . 
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given by 

N = [Tr(XpX^)] ~ l . (54) 

It can be then shown that the generalized decoherence 
function (|52| is normalizable and hermitian, has positive 
diagonal elements, and satisfies the superposition princi- 
ple, provided that [13| 

D(a',a)= 

a'£6c f a£a 

for all coarse grainings {a} of {a}. Thus, this function 
satisfies all consistency tests. 

In the sum-over-histories formulation of gravitational 
systems the usual probability function W is replaced for a 



probability function for a given set of alternatives a, p{a), 
which can be obtained from the decoherence function by 
using the relation D(a',a) rs 5 a > a p(a). We shall take 
these probabilities p{a) as the physical quantities that 
replace quantum states in our non-causal system. So, we 
are here particularly interested in calculating from Eq. 
(15^1) the probability p(a, E; to, a) of a set of alternatives 
for particular values of to, a and P a [B(to,a)], that dis- 
tinguish only the scale factor values on the whole pieces 
of surfaces , B(to,a), which should be spacelike sepa- 
rated from the nonfoliating region. The most general 
expression for the probability p(to, a) (in which we have 
specialized at the particular slicings t = to and a = a) 
that can be obtained from the decoherence function ([5^)) 
is 



p (a, E'; t , a) = p (a, £"; to, a) 

= NTv {XP a [B(t Q , a)} P P a [B{to, a)] X^} = NTr {P a [B(t Q , a)} XpX^P a [B(t , a)] } . (55) 

I 



Finally, the quantity 

f°° f°° / \ 
p(a, E') = p(a, E") = / dto / dap ( a, £^; to, a 1 , 

where the superscript (v) denotes either ' or ", is the most 
general probability for the set of all alternatives that are 
able to distinguish the scale factor values on spacelike 
sections that are spacelike separated from the identified 
nonfoliating region. It will be here regarded to be the 
quantity that replaces the quantum state for a parallel 
phantom universe in Hamiltonian quantum cosmology. 

Because of the non-unitarity of the evolution in the 
neighborhood of the big rip singularity, described in Eq. 
(|55[) by the non-unitary operator X, the probabilities 
given by Eq. (fST)|) for alternative histories completely 
defined on a current local piece located around e.g. our 
galaxy on a given hypersurface, would depend on the 
state of the acronal region [l3j , in the same way as prob- 
abilities for retrodiction histories [1] depend on the cur- 
rent and initial states of the universe. In fact, by the 
cyclic property of the trace, we can rewrite Eq. (|55p as 

p (a, £'; t , a) = iVTr {p f P a [B(t Q , a)} P P a [B(t , a)}} . 

(57) 

where the density matrix, pf, is given by pf — XX' . 
It would mean that, experiments in local laboratories, 
say the solar system or the Galaxy, might give different 
results depending on whether this non-chronal region in 
our future exists or not. So, at least from a theoreti- 
cal point of view, it could glimpse the idea of measuring 
global properties from local experiments. 



V. CONCLUSIONS 

A multiverse scenario can arise in the realm of a phan- 
tom energy dominated universe when we consider the 
complete range of the time interval, smoothing some- 
how the big rip singularity. That multiverse scenario 
comes up from the constraint that we need to impose 
onto the equation of state parameter, in order to obtain 
well-defined values for the scale factor at times after the 
singularity. 

Analytic solutions can be obtained for the quantum 
state of a dark energy dominated universe. If the state 
of the universe is given by a wave function, i.e., if it is a 
pure state, its quantum representation can be expressed 
as a linear combination of Bessel functions and, impos- 
ing the appropriate boundary conditions, we can recover 
the semiclassical approximation in the usual way. Never- 
theless, the most general quantum state for the universe 
should be given by a density matrix, which not with un- 
derstanding suffers from the usual divergence shortcom- 
ings. In order to compute an explicit expression for the 
density matrix in the case being considered, we use a 
particular gauge in which the Hamiltonian cigenfunctions 
can be found, and that can be employed as a basis for 
the space of functions that the Hamiltonian acts upon, 
although this basis is not orthogonal. Then, we found an 
orthonormal set in terms of which the Hamiltonian eigen- 
functions can be expressed and, thus, several expressions 
for the density matrix are given. 

We show that cosmic entangled states between uni- 
verses can take place in the realm of the phantom multi- 
verse. We also give a quantum description of a phantom 
universe when the singularity is replaced for a bounded 
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non-chronal region. In such a case, the generalized quan- 
tum theory is applied and consistent expressions for the 
probabilities of alternative histories are given. 

Although we have succeeded in obtaining a function 
that replaces the conventional notion of quantum state in 
a cosmological spacetime endowed with a bounded multi- 
ply connected region, the model considered in this paper 
is not realistic enough for at least the following reason. 
We have not specifically introduced any matter fields in 
the model, so that this should at best be considered as 
an asymptotic idealization. 

Now some kinds of quantum communication channels 
could be conceived which related the different universes 
that belong to the multiverse. Whether these commu- 
nications could be implemented physically between ad- 
vanced civilizations existing in such universes is a matter 
that requires further consideration. 
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APPENDIX A: ORTHOGONALITY PROPERTIES 
OF BESSEL FUNCTIONS 



The usual formula for the integrals of Bessel functions 
can be obtained from the standard bibliography, to be 0] 



J 



F(z) 



1 



J ll {kt)J v {kt)dt 



fi 2 — v 2 



{kz[J„ +1 (kt)J v (kt) - J^k^Ju+^kt)] - (ji-v)J M {kt)J v (Jst)} (Al) 



r 



Then, a definite integral over t summing from to oo 
can be thought of as the substraction of the two following 
limits, 

,0O J 

/ dt-JJkt)J u (kt) = lim F(z)- limF(z). (A2) 

Jq t z^oo z^O 



Let us compute lim a _vo F(z) first. In this case, taking 
the asymptotic limits for the Bessel's functions, we have 



lim F(z) 

z— »0 



(l 2 — V 2 



{kz 



kz 



2 

fl+V 

1 



zz 

2 J - v){n + v) 

i 



ro + i)r> + 1) n+u 



i 



i 



r(/i + 2)T(u + 1) r(/i + i)i> + 2) 



T{n + l)T(u + 1) 



} 



(ji-v) 



kz 

2) r( M + 2)r(^ + 2) 
(p + v > 0). 



+ (l*-v) 



(A3) 



r 



We can check that Eq. I|A3[) vanishes in the limit for any off at the Planck length is introduced to regularize the 
(/i + v > 0), but it should be taken into account if a cut zero mode Bessel function. For the upper limit we obtain 
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lim F(z) 



/i 2 — V 2 7T 



1 2fc r 1 7T 1 7T 

- cos(z - i/i7r - j) cos(z - ^(i/ + 1)tt - |)] + O(-) 
k sin(§(^ — ju)) 



1 



(A4) 



For + v > 0, it gives: 
1) if /i = ^ 



lim F(^) = 



A* + f 

2) if /j. ^= v and v — \i = 2n (even) 

lim F(z) = 

2 — >00 

3) if ^ v and f — /x = 2n + 1 (odd) 

lim F(z) = — ■ ; 



7T (i/-ji)(i/ + /z) 



(A5) 



(A6) 



(A7) 



Essentially, these are the values of the integral Eq. (|A2[) 
since the limit to zero vanishes. It is left the case for the 
zero mode. We can regularize it by taking some cut-off 
or minimum length, and evaluate the limit 



lim lim 



oo 

dt-J^k^Jvikt), (A8) 



expanding in /i and v the limits Eq. (|A3[) and Eq. (|A4[) . 
i.e. 



1 sin(§^-,i)) 



M + v 



2 



M + v 



{1 + 0{{h-v) 2 )) (A9) 



and 



)1+V 



« -^(l+( M+i /) ln(^)+0(( M+i ,) 2 )) (A10) 

[1 + 1/ fJL + V 2 



so, the limit in Eq. (|A8|) gives 



lim lim 



1 2 
di - J M (M) J„(kt) = lim ln(-). (All) 
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